September 2009 



Evidence for mirror dark matter from 
the CDMS low energy electron recoil spectrum 



R. Foot^ 

School of Physics, 
University of Melbourne, 
Victoria 3010 Australia 



We point out that mirror dark matter predicts low energy ~ 2 keV) electron 
recoils from mirror electron scattering as well as nuclear recoils from mirror ion 
scattering. The former effect is examined and applied to the recently released low 
energy electron recoil data from the CDMS collaboration. We speculate that the 
sharp rise in electron recoils seen in CDMS below 2 keV might be due to mirror 
electron scattering and show that the parameters suggested by the data are roughly 
consistent with the mirror dark matter explanation of the annual modulation signal 
observed in the DAMA/Libra and DAMA/Nal experiments. Thus, the CDMS data 
offer tentative evidence supporting the mirror dark matter explanation of the DAMA 
experiments, which can be more rigorously checked by future low energy electron 
recoil measurements. 
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A successful dark matter theory emerges if one postulates that the fundamental 
laws governing the interactions of the dark matter sector are identical to the ordinary 
matter sector. That is, we assume the existence of a sector which is an exact 
duplicate of the ordinary matter sector, so that the fundamental Lagrangian is: 



Such a dark matter theory can be motivated from simplicity and minimality. It can 
also be motivated from a symmetry reason if left and right handed chiral fields are 
interchanged in the extra sector. This is because the theory then exhibits the space- 
time Z2 parity symmetry, x —x, where chiral left (right) handed ordinary fermions 
transform into chiral right (left) handed mirror fermions; and ordinary bosonic fields 
transform into mirror bosonic fields (see rcf.[l] for the precise transformation). Since 
this theory enlarges the space-time symmetry to include parity (as well as the other 
improper Lorentz transformations), we refer to the particles in the extra sector as 
mirror particles. The standard model extended with a mirror sector was first studied 
in ref. [1] and shown to be a phenomenologically consistent renormalizable theory (for 
a review and more complete list of references see ref. [2]) ^. 

If we include all interaction terms consistent with renormalizability and the sym- 
metries of the theory then we must add to the Lagrangian a U{1) kinetic mixing 
interaction [4] and Higgs - mirror Higgs quartic coupling[l]: 



where F^i, {F'^i) is the ordinary (mirror) ^7(1) gauge boson field strength tensor 
and ((/)') is the electroweak Higgs (mirror Higgs) field. The most general Higgs 
potential, including the quartic Higgs mixing term (above) was studied in ref. [1] and 
shown to have the vacuum (0) = (0') for a large range of parameters. With this 
vacuum, the masses of the mirror particles are all identical to their ordinary matter 
counterparts. Note that interactions between ordinary and mirror particles at low 
energies are expected to be dominated by the U{1) kinetic mixing interaction. This 
interaction leads to photon - mirror photon kinetic mixing, which provides the key 
to testing this theory (for a review, see e.g. ref. [2]). 

In this framework, dark matter is comprised of stable massive mirror particles: 
e', if'. He' , O' ... etc, with known masses. Dark matter in the galactic halo is then 
a spherically distributed self interacting mirror particle plasma which can be stable 
from collapse provided that a heat source(s) exists. In fact, ordinary supernova can 
plausibly supply the required heating if photon-mirror photon are kinctically mixed 
with e ~ 10~^[5] ^. For kinetic mixing of this magnitude about half of the total 

^Note that successful big bang nucleosynthesis (BBN) and large scale structure (LSS) requires 
effectively asymmetric initial conditions in the early Universe, T' <C T" and nt'/nb « 5. See ref. [3] 
for further discussions. 

mirror sector with such kinetic mixing is consistent with all known laboratory, astrophysical 
and cosmological constraints [6, 7, 8, 9, 10] . 



^ = CsM{e, u, d, 7, ...) + CsM{e', u', d', 7', ...) 



(1) 



(2) 
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energy emitted in ordinary Type II Supernova explosions (~ 3 x 10 erg) will be 
in the form of light mirror particles {u'^^^, e'^,'~f') originating from kinetic mixing 
induced plasmon decay into e'~^e'~ in the supernova core[10]. This implies a heating 
of the halo (principally due to the e'^ component), of around: 

Lilt-in ~ ^ X 3 X 10^^ er^^^^i^ ~ 10^^ erg/s for Milky Way (3) 

It turns out that this matches (to within uncertainties) the energy lost from the halo 
due to radiative cooling[5]: 

L'::^ergy-out = A /^^ ulAnr'dr - 10^^ erg/s for Milky Way. (4) 

In other words, a gaseous mirror particle halo can potentially survive without col- 
lapsing because the energy lost due to dissipative interactions can be replaced by 
the energy from ordinary supernova explosions. Presumably there is some detailed 
dynamical reasons maintaining this balance, which of course, may be difficult to 
elucidate due to the complexity of that particular problem. 

Kinetic mixing of around e ~ 10"^ is also implicatcd[ll] by the annual modula- 
tion signal found in the DAMA/Nal and DAMA/Libra experiments[12]. In ref.[ll] 
(updating and improving earlier studies[13]), it was shown that elastic scattering of 
the heavy mirror oxygen {O') component off the target nuclei can nicely explain the 
DAMA signal if e ~ 10~^ including the observed recoil energy dependence of the 
annual modulation amplitude. Importantly, this explanation is also consistent with 
the null results of the other direct detection experiments. 

In view of the success of the mirror dark matter theory in explaining the existing 
direct detection experiments, it is important to look for other ways to further test 
this theory. In this note we would like to point out another signature of mirror 
dark matter relevant for direct dark matter detection experiments, which has not 
been previously discussed. In addition to target nuclear recoils, low energy electron 
recoils are also predicted arising from mirror electron ordinary electron scattering. In 
order to compute the mirror electron interaction rate, we need to know the number 
density of mirror electron's at the Earth's location, their velocity distribution, and 
interaction cross section. 

Assuming that the dark matter density is^ u = 0.3 GeV/cm^ at the earth's 
location, and that the halo mass is dominated by the H', He' component [with mass 
fraction Yhc' = nHe'i^He/ {'n'H''fnH + ^T-iJe'^^JJe)]; we expect a mirror electron number 
density, ^e', of: 

„, = ^fl-%) (6) 



rrip 



where rup is the proton mass. 

*We use natural units h = c=l unless otherwise stated. 
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The halo mirror particles arc i)iesumcd to form a self interacting spherically 
distributed plasma at temperature T. [This is necessary to explain the flat rotation 
curves in spiral galaxies]. The dynamics of the mirror particle plasma has been 
investigated previously[5, 13], where it was found that the condition of hydrostatic 
equilibrium implied that the temperature of the plasma satisfied: 

T = \mvlt , (6) 

where in — X^rijmi/ X^n, [i — e', H', He', O'....] is the mean mass of the particles in 
the plasma, and Vj-ot ~ 254 km/s is the rotational velocity of the sun around the 
center of the galaxy[14]. Assuming the plasma is completely ionized, a reasonable 
approximation since in turns out that the temperature of the plasma is ~ | keV 
(see below) we find: 

^ = I (7) 

f^p 2 — jYhe' 

Hence from Eq.(6) we find the temperature to be: 

T 0.48 keV for Yne' = 1 

T 0.18 keV for Yne' = (8) 

Thus the temperature is quite sensitive to Yhc', but mirror BBN calculations sug- 
gest Yffe' ~ 0.9 for e ~ 10~^[15]. Anyway, the gas of mirror particles will have a 
Maxwellian velocity distribution. 



fiiv) 
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„2 /„2 



= e-^ /"ol^l (9) 



where the index % labels the particle type [i = e', if', ife', O', ...]. 

Clearly, the velocity dispersion of the particles in the mirror matter halo depends 
on the particular particle species and from Eq.(6) satisfies: 

Vl^\ = V'rot- (10) 

rrii 

Note that if m,j ^ m, then v^li] <^ v^^^. Consequently heavier mirror particles have 
their velocities (and hence energies) relative to the earth boosted by the Earth's 
(mean) rotational velocity around the galactic center, Vrot- This allows a heavy 
mirror particle in the mass range ~ mo ~ 15 GeV to provide a significant annual 
modulation signal in the energy region probed by DAMA (2 < E^/keV < 6), which 
it turns out has the right properties to fully account [11] for the data presented by 
the DAMA collaboration[12]. On the other hand the mirror electron, being much 
lighter than the mean mass, m, will have a very large velocity dispersion: 

4777' 

'^oi^') — '^lotT. — - for He' dominated halo, i.e. Yne' = 1 (H) 

oTflp 
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That is, we expect vo{e') ~ 12, 000 km/s for Yhs' ~ 1- Note that at these velocities, 
we can, to a good approximation, neglect the motion of the Earth through the halo 
as far as the e' component is concerned. Thus, since T ~ \keV , we expect an 

electron scattering signal to be confined to low energy electron recoils ~ 2 keV. 

We now discuss the interaction cross section between ordinary and mirror elec- 
trons. The effect of the photon-mirror photon kinetic mixing is to induce a small 
coupling of ordinary photons to mirror electrons, of magnitude ee. This enables 
mirror electrons to Rutherford scatter off free ordinary electrons with cross section 
in the non-relativistic limit: 

^'^ ^ (12) 



where 

A ^ ^ (13) 

Here Er is the recoil energy of the target electron, initially presumed at rest relative 
to the incoming mirror electron of velocity v. This cross section should also approx- 
imate the inelastic scattering of halo mirror electrons off bound atomic electrons 
provided that the recoil energies are much larger than the binding energy. 

Recently, the CDMS collaboration has released data on low energy electron scat- 
tering down to a recoil energy of 1 keV on a Germanium target [17, 18], which we 
reproduce in figure 1. Such an experiment is potentially sensitive to mirror electron- 
ordinary electron scattering, where we expect a sharp rise in interaction rate below 
2 keV. Interestingly, such a rise in interaction rate was seen in that experiment, 
although it should also be noted that this rise is near the CDMS energy thresh- 
old. Certainly, any conclusions drawn from an analysis of the Eji < 2 keV data 
are tentative and will need to be confirmed by future low threshold electron recoil 
measurements. 
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Figure 1: Efficiency corrected low-energy electron-scattering spectrum as measured 
by CDMS/Ge[17]. The dashed line is a simple background model. 



In principle, it would be desirable to properly account for the cross section of 
mirror electrons on bound atomic electrons which is a non-trivial quantum mechan- 
ical problem and is beyond the scope of this initial study. In this paper, we will 
perform a crude analysis using the following simple approximation. We consider the 
scattering only on the loosely bound Ge electrons in the 3d, 4s and 4p outer shells, 
which have atomic binding energy ~ 30 eV[16]. There are 14 such electrons per Ge 
atom. The remaining 18 electrons occupy the more tightly bound inner shells all of 
which have binding energy ~ 120 eV (and the 10 inner most electrons have binding 
energy ~ 1.2 keV) and have typically much faster velocities, comparable to Vo{e'). 
We expect the interactions on these inner shell electrons to be suppressed relative to 
the more weakly bound outer shell electrons. Furthermore, we approximate the 14 
weakly bound outer shell electrons of Ge, as free and at rest. We expect the result- 
ing cross section to be valid to within around 20 — 30% for the recoil energy range 
1 < Eji/keV < 2. With this approximation, the predicted differential interaction 
rate is: 

dR f da fe'{v) , I 3 
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- gNTUe'-^ / -Tryd'v (14) 

(Er) k\v\ 

where Nt is the number of target atoms per kg of detector and k — [7rvQ(e')]^/^ is the 
MaxweUian distribution normahzation factor. The quantity g = 14, is the number 
of loosely bound atomic electrons in Ge as we discussed above. Note that the lower 
velocity limit, VminiEn), is given by the kinematic relation: 



R 



V rUe 

The velocity integral in Eq.(14) can be analytically solved leading to: 

dR A / 2e-^" 



(15) 



where x = Vmin/vo- 

In order to compare with the experimentally measured rate, we must convolve 
this rate, with a Gaussian to take into account the finite detector resolution: 

dR 1 r dR 



dE^ 



= ^ / (17) 
(Tv27r J dEji 



Here E]l is the 'measured recoil energy' while actual recoil energy we denote as Eji. 
The detector averaged resolution, u, is measured to be[18]: 

cr = ^(0.293)2 + {0.056y Er keV . (18) 

We also incorporate a detection efficiency effect as follows. Prom figure 2 of ref.[18] 
it is evident that the detection efficiency goes to zero at Er 0.8 keV, and this 
suggests a lower recoil energy limit of around 0.8 keV. That is, we integrate Eq.(17) 
down to actual recoils of 0.8 keV. Aside from this, it is unnecessary to include 
the detection efficiency explicitly in Eq.(17) since the data presented is efficiency 
corrected. 

Numerically, we can now compare the predicted rate described above with the 
low energy recoil electron scattering data from the CDMS/Ge experiment. The rate 
depends on the parameters e,YHe'- We fix Yhc' ~ 0.9 which is the value predicted 
from mirror BBN analysis[15]. In figure 2 we give our results for e = 6 x lO"^'^ 
(dashed line), e = 7 x 10"^° (solid fine) and e = 8 x 10"^° (dotted line). We have 
assumed a nearly flat background rate of R(background) = 1.9 — 0.09 x E^ (shown 
as the dashed line in figure 1), which is similar to the approach of ref.[17]. 
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Figure 2: Mirror dark matter induced electron recoils + background model, for the 
parameters Yhc' = 0.9 and e = 6 x 10~^° (dashed line), e = 7 x 10~^° (solid line) 
and e = 8 X 10~^° (dotted line), compared with the CDMS/Ge data. 

Anyway, despite the various systematic uncertainties and other deficiencies in 
our analysis, it is nevertheless quite interesting that the value of e suggested by the 
data is around e 7 x 10~^°. This value is consistent with the mirror dark matter 
explanation of the DAMA annual modulation signal which requires: 



where ^o' = no'Tno/cu is the proportion of O' by mass in the galactic halo. If taken 
seriously the CDMS data suggest a relatively large mirror metal {O') proportion in 
the halo. Such a scenario is possible given the asymmetric evolution of ordinary 
and mirror matter implied by the asymmetric initial condition T' <^T inferred to 
exist in the early Universe (sec e.g. ref. [3] for discussions). For example T' <^ T 
in the early Universe implies a higher Yhc' fraction (c.f Y^e) which leads to more 
rapid mirror star evolution and early structure formation in the mirror sector[19]. 
The more rapid stellar evolution in the mirror sector might be expected to lead to 
a higher metal proportion c.f. the ordinary matter sector. 

Finally in figure 3 we examine the effect of changing the Yhc' parameter for fixed 
e. Fixing e = 7 x 10~^°, we show results for Yhc' — 0.3, 0.6, 0.9. As the figure 




(19) 
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demonstrates the results dependent somewhat sensitively on the value of Yhc' which 
is due to the Temperature dependence (or equivalently, i'o(e') dependence) on Yhc' 
as discussed earlier. 
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Figure 3: Mirror dark matter induced electron recoils + background model, for fixed 
e = 7 X 10-1° and Yhc' = 0.3 (dotted line), Y^e' = 0.6 (dashed line) and Yhc' = 0.9 
(solid fine). 

Recall that the mirror dark matter interpretation of the DAMA annual modula- 
tion signal depends less sensitively on l//e'[ll]- Thus, ultimately low recoil energy 
electron scattering measurements might provide a useful way to measure Yhc', which 
can then be checked against the theoretical expectation from mirror BBN analysis, 
which suggest Y^e' ~ 0.9[15]. 

In conclusion, we have pointed out that mirror dark matter predicts low energy 
electron recoils (Er ^ 2 keV) from mirror electron scattering as well as nuclear 
recoils from mirror ion scattering. The former effect is examined and applied to the 
recently released low energy electron recoil data from the CDMS collaboration. We 
speculate that the sharp rise in electron recoils seen in CDMS below 2 kcV might 
be due to mirror electron scattering and show that the parameters suggested by the 
data are roughly consistent with the mirror dark matter explanation of the annual 
modulation signal observed in the DAMA/Libra and DAMA/Nal experiments. This 
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interpretation of the CDMS data can be more rigorously checked by future low 
energy electron recoil measurements, which we await with interest. 
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